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Abstract. 

We review how our current understanding of the light element synthesis during 
the Big Bang Nucleosynthesis era may help shed light on the identity of particle dark 
matter. 



1. Introduction 

In the late 40s and throughout the 50s a number of visionary scientists including Alpher, 
Fermi, Follin, Gamow, Hayashi, Herman, and Turkevich attempted to explain nuclear 
abundance patterns observed in the nearby Universe, such as the peculiar high helium 
mass fraction Y p « 0.25. This initially speculative work on an era of nucleosynthesis 
(element formation) in an expanding Universe at very high temperature T ~ 10 9 K 
developed slowly but steadily over the coming decades into what is now known as the 
standard model of Big Bang nucleosynthesis (BBN). The idea that the Universe may 
have undergone a very hot and dense early phase got triggered by the observations 
of Hubble in the 1920s, of the recession velocity of galaxies being proportional to their 
inferred distance from the Milky Way, which were most elegantly explained by a Universe 
in expansion. The "expanding, hot Big Bang" idea received further support by the 
observation of the cosmic microwave background radiation (CMBR) by Penzias and 
Wilson in 1965, believed to be the left-over radiation of the early Universe. Detailed 
observational and theoretical studies of BBN as well as the CMBR and the Hubble flow 
have developed into the main pillars on which present day cosmology rests. 

BBN takes place between eras with (CMBR) temperatures T ~ 3 MeV and 
T ~ lOkeV, in the cosmic time window t ~ 0.1 — 10 4 sec, and may be characterized as 
a freeze-out from nuclear statistical equilibrium of a cosmic plasma at very low ~ 10~ 9 
baryon-to-photon number ratio (cf. Section [2]), conditions which are not encountered 
in stars. It produces the bulk of 4 He and 2 H (D), as well as good fractions of 3 He 
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and 7 Li observed in the current Universe, whereas all other elements are believed 
to be produced either by stars or cosmic rays. In its standard version it assumes a 
Universe expanding according the laws of general relativity, at a given homogeneously 
distributed baryon-to-photon ratio r] b , with only Standard Model particle degrees of 
freedom excited, with negligible /j| <T lepton chemical potentials, and in the absence 
of any significant perturbations from primordial black holes, decaying particles, etc. 
By a detailed comparison of observationally inferred abundances (cf. Section [3]) with 
those theoretically predicted, fairly precise constraints/conclusions about the cosmic 
conditions during the BBN era may thus be derived. BBN has been instrumental, 
for example, in constraining the contribution of extra "degrees of freedom" excited 
in the early Universe to the total energy density, such as predicted in many models 
of particle physics beyond the standard model. Such contributions may lead to an 
enhanced expansion rate at T ~ 1 MeV implying an increased 4 He mass fraction. It 
is now known that aside from baryons and other subdominant components not much 
more than the already known relativistic degrees of freedom (i.e. photons 7's, electrons 
and positrons e ±, s, and three left-handed neutrinos z/s) could have been present during 
the BBN era. BBN is also capable of constraining very sensitively any non-thermal 
perturbations as induced, for example, by the residual annihilation of weak scale dark 
matter particles (Section EJ), or by the decay of relic particles (Section [7|) and the possible 
concomitant production of dark matter. Moreover, the sheer presence of negatively 
charged or strongly-interacting weak mass-scale particles during BBN (Section [6]) may 
lead to dramatic shifts in yields of light element through the catalytic phenomena. BBN 
may therefore constrain properties and production mechanisms of dark matter particles, 
and this chapter aims at revealing this connection. 

It is possible that the biggest contribution of BBN towards understanding the 
dark matter enigma has already been made. Before the advent of precise estimates 
of the fractional contribution of baryons to the present critical density, Qb ~ (0.02273 ± 
0.00062)//i 2 , where h is the Hubble constant in units 100 km s _1 Mpc _1 , by detailed 
observations and interpretations of the anisotropics in the CMBR [I], BBN was the 
only comparatively precise mean to estimate Q^. As it was not clear if the "missing" 
dark matter was simply in form of brown dwarfs, white dwarfs, black holes (formed 
from baryons), and/or T ~ 10 6 K hot gas, various attempts to reconcile a BBN era at 
large ~ 1 with the observationally inferred light element abundances were made. 
These included, for example, BBN in a baryon-inhomogeneous environment, left over 
possibly due to a first-order QCD phase transition at T m 100 MeV, or BBN with 
late-decaying particles, such as the supersymmetric gravitino (for reviews, cf. [21 [HE]). 
Only continuous theoretical efforts of this sort, and their constant "failure" to account 
for large f^, gave way to the notion that the dark matter must be in form of "exotic", 
non-baryonic material, such as a new fundamental particle investigated in the present 
book. 
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2. Standard BBN - theory 

Standard BBN (SBBN) theory is well understood and described in detail in many 
modern cosmology text books. The essence of SBBN is represented by a set of 
Boltzmann equations that may be written in the following schematic form: 

^ = ~ H ( T ) T ^§ = E( r ^ + W« + •••), (i) 

where Y% = n%/s are the time t (or temperature T) dependent ratios between number 
density and entropy density s of light elements % = 1 H, n, D, 4 He, etc.; the r^.., 
are generalized rates for element interconversion and decay that can be estimated by 
experiments and/or theoretical calculations, and H(T) is the temperature-dependent 
Hubble expansion rate. The system of equations (JTJ) assumes thermal equilibrium 
e.g. Maxwell-Boltzmann distributions for nuclei, which is an excellent approximation 
maintained by frequent interactions with the numerous 7's and e ±, s in the plasma. The 
initial conditions for this set of equations are well-specified: for temperatures much in 
excess of the neutron-proton mass difference, neutron and proton abundances are equal 
and related to the baryon to entropy ratio, Y neutron ~ Y prot0 n — ^ n baryon/s, while the 
abundance of all other elements is essentially zero. At temperatures relevant to BBN, 
the baryonic contribution to the Hubble rate is minuscule, and H(T) is given by the 
standard radiation-domination formula: 

ffm ^2 v ( 8n 3 g*G N \ 1/2 7 

H{T) = T X I — I , where g* = g boson + -gf erm ion , (2) 

with the g's denoting the excited relativistic degrees of freedom. This expression 
needs to be interpolated in a known way across the brief epoch of the electron- 
positron annihilation, after which the photons become slightly hotter than neutrinos 
and H(T) ~ T 9 2 /(178 s), where Tg is the photon temperature in units of 10 9 K. A number 
of well-developed integration routine that go back to an important work of Wagoner, 
Fowler, and Hoyle [5], allow to solve the BBN system of equations numerically and 
obtain the freeze-out values of the light elements. A qualitative " computer-free" insight 
to these solutions can be found in e.g. Ref. [6]. 

In a nutshell, SBBN may be described as follows. After all weak rates fall below 
the Hubble expansion rate, the neutron-to-proton ratio freezes to ~ 1/6, subject to a 
slow further decrease to ~ 1/7 by Tg ~ 0.85 via neutron decay and out-of-equilibrium 
weak conversion. At this point, to a good approximation, all neutrons available will be 
incorporated into 4 He, since it is the light element with the highest binding energy per 
nucleon. Synthesis of 4 He, and all other elements, has to await the presence of significant 
amounts of D (the "deuterium bottleneck"). This occurs rather late, at Tg ~ 0.85, since 
at higher Tg the fragile D is rapidly photodisintegrated by the multitude of CMBR 
photons. At Tg < 0.85 the fairly complete nuclear burning of all D then results in only 
trace amounts O(10 -5 ) of D (and 3 He) being leftover after SBBN has ended. Elements 
with nucleon number A > 4 are even less produced due to appreciable Coulomb barrier 
suppression at such low Tg, resulting in only O(10~ 10 ) of 7 Li, and abundances of other 
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isotopes even lower. SBBN terminates due to the combination of a lack of free neutrons 
and the importance of Coulomb barriers at low T. In the following a few more details 
are given: 

O(0. 1) abundances: 4 He. The 4 He mass fraction Y p is dependent on the timing of 
major BBN events, such as the neutron-to-proton freeze-out at T ~ 0.7 MeV, post- 
freezeout neutron depletion before the deuterium bottleneck, and the position of this 
bottleneck itself as a function of temperature. Consequently, Y p is dependent on such 
well-measured quantities as Newton's constant, the neutron-proton mass difference, 
neutron lifetime, deuterium binding energy, and to a much lesser degree on less precisely 
known values for the nuclear reaction rates. This sensitivity to the timing of the BBN 
events makes 4 He an important probe of the Hubble expansion rate, and of all possible 
additional non-standard contributions that could modify it. The SBBN predicts Y p with 
an impressive precision, Y p = 0.2486 ± 0.0002, where we use the most recent evaluation 

® 

O(10 -5 ) abundances: D and 3 He. Deuterium and 3 He BBN predictions are 
more sensitive both to nuclear physics and to rjt, input. Reactions involving these 
elements are well measured, and with the current WMAP input SBBN is capable of 
making fairly precise predictions of these abundances: D/H = 2.49 ± 0.17 x 10~ 5 ; 
3 He/H = (1.00 ± 0.07) x 10~ 5 . 

O(10 -10 ) abundances: 7 Li. Among all observable BBN abundances, 7 Li is the most 
sensitive to the i] b and nuclear physics inputs. The actual observable that BBN predicts 
is the combined abundance of 7 Li and 7 Be, as later in the course of the cosmological 
evolution 7 Be is transformed into 7 Li via electron capture. At the CMBR-measured value 
of the baryon-to-photon ratio rjb, more than 90% of primordial lithium is produced in 
the form of 7 Be in the radiative capture process, 4 He+ 3 He— > 7 Be+7. As the rate for this 
process per each 3 He nucleus is much slower than the Hubble rate, the output of 7 Be is 
almost linearly dependent on the corresponding S"-factor for this reaction. With recent 
improvement in its experimental determination |8j [9j [10], the current ~15% accuracy 
prediction for 7 Be+ 7 Li stands at 5.2At° H x 10" 10 [7j. 

O(10 -14 ) and less abundances: 6 Li and A > 9 elements. 6 Li is formed in the BBN 
reaction 

4 He + D -> 6 Li + 7, Q = 1.47MeV (3) 

which at BBN temperatures is ~ four orders of magnitude suppressed relative to other 
radiative capture reactions such as 4 He+ 3 H-^ 7 Li + 7, and ~ seven/eight orders of 
magnitude suppressed relative to other photonless nuclear rates. The reason for the 
extra suppression is in a way accidental: it comes from the same charge to mass ratio for 
4 He and D, which inhibits the El transition, making this radiative capture extremely 
inefficient. This results in O(10~ 14 ) level prediction for primordial 6 Li which is well 
below the detection capabilities. Heavier elements with A > 9 such as 9 Be, 10 B and 
11 B are never made in any significant quantities in the SBBN framework, and the main 
reason for that is the absence of stable A = 8 nuclei, as 8 Be is underbound by 92 keV 
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and decays to two a. 

3. Observed light element abundances 

In the following, we will briefly discuss the observationally inferred light elements, 
element by element. Here the discussion will also include isotopes, or isotope ratios, 
such as 6 Li, 9 Be, and 3 He/ 2 H, which are not always considered in SBBN, but which are 
very useful to constrain deviations from SBBN. 

3.1. 4 He 

The primordial 4 He/H ratio is inferred from observations of hydrogen- and helium- 
emission lines in extragalactic low-metallicity Hll-regions and compact blue galaxies, 
illuminated by young star clusters. Two particular groups have performed such 
analysis for years now, with their most recent results Y p « 0.2477 ± 0.0029 |TT] and 
Y p m 0.2516 ± 0.0011 [12]. These estimates are (surprisingly) considerably larger than 
earlier estimates by both groups (i.e. 0.239 and 0.242, respectively), explained in large 
parts by a new estimate for Hel emissivities [13]. Other differences with respect to 
older studies, and/or between the two new studies themselves, are the adopted rates for 
collisional excitation of H- (He-) emission lines, corrections for a temperature structure 
in these galaxies ("temperature variations"), corrections for the presence of neutral 4 He 
("icf - ionisation corrections"), as well as corrections for throughs in the stellar spectra at 
the position of the 4 He- (H-) emission lines ("underlying stellar absorption"). All of these 
may have impact on the ^ 1% level. This, as well as the comparatively large change from 
earlier estimates (coincidentally going into the direction of agreement with the SBBN 
prediction of Y p « 0.248), implies that a conservative estimate Y p ~ 0.249 ± 0.009 [H] 
(see also Y p m 0.250 ± 0.004) [15] of the error bar, maybe more appropriate when 
constraining perturbations of SBBN. 

3.2. D 

For the observational determination of primordial D/H-ratios high- resolution 
observations of low-metallicity quasar absorption line systems (QALS) are employed 
(cf. [16j [T7J [HJ [SI [20]). QALS are clouds of partially neutral gas which fall on the 
line of sight between the observer and a high-redshift quasar. The neutral component 
in these clouds yields absorption features, for example, at the redshifted position of the 
Lyman-a wavelength. For the very rare QALS of sufficiently simple velocity structure, 
one may compare the absorption at the Lyman-a position of H with that of D (shifted 
by 81 km s _1 ) to infer a D/H ratio. Here the low metallicity of these QALS is conducive 
to make one believe that stellar D destruction in such clouds is negligible. Currently 
there exist only about 6 — 8 QALS with D/H determinations. When averaged they 
yield typical 2.68 - 2.82 x 10" 5 pH EH E21 ED] for the central value, with inferred 
statistical la error bars of 0.2 — 0.3 x 10~ 5 , comparing favorably to the SBBN prediction 
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of 2.49 ± 0.17 x 10~ 5 [7J at the WMAP inferred r]},. Nevertheless, as an important 
cautionary remark, the various inferred D/H ratios in QALS show a spread considerably 
larger than that expected from the above quoted error bars only. This is usually a 
sign of the existence of unknown systematic errors. Until these systematics are better 
understood should primordial values as high as D/H^ 4 x 10~ 5 therefore not considered 
to be ruled out. 



3.3. 3 He/D 

Observational determinations of 3 He/H-ratios are possible within our galaxy which is 
chemically evolved. The chemical evolution of 3 He is, however, rather involved, with 
3 He known to be produced in some stars and destroyed in others. Furthermore, any 
D entering stars will be converted to 3 He by proton burning. The net effect of all this 
is an observed approximate constancy of (D+ 3 He)/H^ 3.6 ± 0.5 x 10 -5 [23] over the 
last few billion years in our galaxy. Whereas the relation of galactic observed 3 He/H 
ratios to the primordial one is obscure, the ratio of 3 He/D as observed in the presolar 
nebulae is invaluable in constraining perturbations of SBBN. This ratio 0. 831^25 [23] 
(where the error bars are obtained when using the independent 2a ranges of 3 He/H and 
D/H) provides a firm upper limit on the primordial 3 He/D [24]. This is because 3 He 
may be either produced or destroyed in stars, while D is always destroyed, such that 
the cosmological 3 He/D ratio may only grow in time. 



3.4. 7 Li 

7 Li/H ratios may be inferred from observations of absorption lines (such as the 6708A 
doublet) in the atmospheres of low-metallicity galactic halo stars. When this is done 
for stars at low metallicity [Z], 7 Li/H ratios show a well-known anomaly (with respect 
to other elements), i.e. 7 Li/H ratios are constant over a wide range of (low) [Z] and 
some range of temperature (the "Spite plateau"). As most elements are produced by 
stars and/or cosmic rays, which themselves produce metallicity, the 7 Li Spite plateau is 
believed to be an indication of a primordial origin of this isotope. This interpretation is 
strengthened by the absence of any observed scatter in the 7 Li abundance for such 
stars. There have been several observational determinations of the 7 Li abundance 
on the Spite plateau. Most of them fall in the range 7 Li/H^ 1 — 2 x 10~ 10 such as 
1.23l{j:32 x 10~ 10 [251 [26] and 1.1 - 1.5 x 10~ 10 [27], with some being somewhat higher 
such as 2.19 ± 0.28 x 10~ 10 [28]. Here differences may be due to differing methods of 
atmospheric temperature estimation. These values should be compared to the SBBN 
prediction 5.24lg;g7 x 10~ 10 [7] (with la error estimates), clearly indicating a conflict 
which is often referred to as the " lithium problem" . It is essentially ruled out that this 
problem be solved by, either, an erroneous atmospheric temperature determination, or 
significant changes in 7 Li producing/destroying nuclear SBBN rates. There remain only 
two viable possibilities of a resolution to this statistically significant (4 — 5)a problem. 
First, it is conceivable that atmospheric 7 Li has been partially destroyed in such stars due 
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to nuclear burning in the stellar interior. Though far from understood, one may indeed 
construct (currently ad hoc) models which deplete 7 Li by a factor « 2 in such stars, 
while respecting all other observations [29j EQ]. Second, it is possible that the lithium 
problem points directly towards physics beyond the SBBN model, possibly connected 
to the production of the dark matter (cf. Section [7J. 

3.5. & Li and 9 Be 

The isotope of 6 Li is usually not associated with BBN, as its standard BBN production 
6 Li/H~ 10~ 14 is very low. However, the smallest deviations from SBBN may already lead 
to important cosmological 6 Li abundances. It is therefore interesting that the existence 
of 6 Li has been claimed in about ~ 10 low-metallicity stars [27], with, nevertheless, 
each of these observations only at the 2 — 4a statistical significance level. Asplund 
et al. infer an average of 6 Li/ 7 Li^ 0.044 (corresponding to 6 Li/H^ 6 x 10~ 12 ) for 
their star sample, whereas Cayrel et al. infer 6 Li/ 7 Li^ 0.052 ± 0.019. for the star 
HD84937 [31]. Such claims, if true, would be of great interest, as the inferred 6 Li/ 7 Li in 
very low metallicity stars is exceedingly hard to explain by cosmic ray production |32j, 
though in situ production in stellar flares may be conceivable [33]. Moreover, the 6 Li 
observations seem to be consistent with a plateau structure at low metallicity as expected 
when originating right from BBN. However, recent work [31] has cast significant shadow 
over the claim of elevated 6 Li/ 7 Li ratios at low [Z]. Similar to 7 Li, 6 Li is inferred from 
observations of atmospheric stellar absorption features. Unlike in the case of D and H 
in QALS, the absorption lines of 7 Li and 6 Li are always blended together. 6 Li/ 7 Li ratios 
may therefore be obtained only by observations of a minute asymmetry in the 6708 line. 
Such asymmetries could be due to 6 Li, but may also be due to asymmetric convective 
motions in the stellar atmospheres. The analysis in Cayrel et al. [31] prefers the latter 
explanation. 

Unlike the case of 6 Li, detection of 9 Be in many stars at low metallicities is not 
controversial. Observations of 9 Be [35, 36J are far above the O(10~ 18 ) SBBN prediction, 
and exhibit linear correlation with oxygen, clearly indicating its secondary (spallation) 
origin [37]. The lowest level of detected 9 Be/H is at ~fewxl0~ 14 , which translates into 
the limit on the primordial fraction at 2 x 10~ 13 [38J, assuming no significant depletion 
of 9 Be in stellar atmospheres. 

4. Cascade nucleosynthesis from energy injection 

The possibility that BBN may be significantly perturbed by the presence of energetic, 
non-thermal SM particles in the plasma has first received detailed attention in the 
1980s [^Ql SH S21 S3J S3, SHI SZ1 - Though much of the pioneering work had been done, 
only recently the first fully realistic calculations of coupled thermal nuclear reactions 
and non-energetic phenomena have been presented [HJ HHJ [50j EE] . Energetic particles 
may be injected as products of the decay or annihilation of relic non-SM particles, or 
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Figure 1. Constraints on the abundance Q,xh 2 of relic particles decaying at tx 
assuming Mx = 100 GeV for the particle mass. The most stringent limits are given, 
from early to late times, by 4 He, D, 6 Li, and 3 He/D overproduction, respectively. The 
various lines are for different logio-Bh, as labeled, where is the hadronic branching 
ratio. ^,From Ref. [39] , 

via perhaps more exotic mechanisms such as evaporation of primordial black holes or 
super symmetric Q-balls. The injected energetic photons 7's, electron/positrons e s, 
neutrinos u's muons /r^s, pions 7r's, nucleons and antinucleons ATs and N's, gauge 
bosons Z's and W s, etc. may be considered as the "cosmic rays" of the early Universe. 
In contrast to their present day counterparts, and with the exception of neutrinos, these 
early cosmic rays thermalize rapidly within a small fraction of the Hubble time H~ l (T) 
for all cosmic temperatures above T ~ 1 eV. This, of course, happens only after all 
unstable species (i.e. 7r's, /i's Z's and W ±, s) have decayed leaving only 7's, e^s, z/'s and 
TV's. Many of the changes in BBN light-element production occur during the course of 
this thermalization. One often distinguishes between hadronically (7r's, iV's, and iV's) 
and electromagnetically (7's, e s) interacting particles, mainly because the former may 
change BBN yields at times as early as r> 0.1 sec (i.e. T<3MeV), whereas the latter 
only have impact for r > 10 5 sec (i.e. T < 3 keV). In the following we summarize the most 
important interactions and outline the impact of such particles on BBN. For hadronically 
interacting particles these effects include: 

(i) n s may cause charge exchange, i.e. it" +p — > tt° + n between 1 MeV > T > 300 keV 
thereby creating extra neutrons after n/p freezeout and increasing the helium mass 
fraction Y p . 

(ii) Antinucleons injected in the primordial plasma preferentially annihilate on 
protons, thereby raising the effective n/p-ratio and increasing Y p . 

(iii) At higher temperatures, neutrons n's completely thermalize through magnetic 
moment scattering on (T ^ 80 keV), whereas protons p's do so through Coulomb 
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interactions with and Thomson scattering off CMBR photons (T> 20keV). Any 
extra neutrons at T ~ 40 keV may lead to an important depletion of 7 Be. 

(iv) At lower temperatures, both, energetic neutrons and protons may spall 4 He, e.g. 
n + 4 He — > 3 H + p + n + (vr's), or n + 4 He — > D +p + In + (vr's). Both reactions are 
important as they may either increase the 2 H abundance or lead to 6 Li formation 
via the secondary non-thermal reactions of energetic 3 H( 3 He) on ambient a's. 

The main features of electromagnetic injection are: 

(i) Energetic 7's may pair-produce on CMBR photons, i.e. 7+7cmbr — ► e~+e + as long 
as their energy is above the threshold Eq ~ m 2 e /22T for this process. The created 
energetic e^ in turn inverse Compton scatter, i.e. e ± + 7cMBR e ± + 7, to produce 
further 7's. Interactions with CMBR photons completely dominate interactions 
with matter due to the exceedingly small cosmic baryon-to-photon ratio rj. 

(ii) Only when E^<Ec do interactions with matter become important. These include 
Bethe-Heitler pair production 7 + p( 4 He) — > p( 4 He) + e + + and Compton 
scattering 7 + e~ — > 7 + e~ off plasma electrons, as well as photodisintegration 
(see below). 

(iii) A small fraction of 7's with E 1 < E c may photodisintegrate first D at T < 3 keV, 
when E c becomes larger that Ef 2.2 MeV, the D binding energy, and later 4 He, at 
T<0.3keV since E 4Hc ks 19.8 MeV. Such processes may cause first, D destruction, 
and later, D and 3 He production, and more importantly 3 He/D overproduction. 
They may also lead to 6 Li production. 

In the context of non-thermal energy injection related to particle dark matter, there 
are two very important processes that have profound impact on 6 Li and 7 Li abundances 
and deserve further comments. Energetic 3 He and 3 H produced via electromagnetic 
or hadronic energy injection (i.e. via spallation or photodisintegration) provide the 
possibility of efficient production of 6 Li via the non-thermal nuclear reactions on thermal 
4 He: 

3 H+ 4 He -> 6 Li+n, Q = -4.78MeV; 3 He+ 4 He -> 6 Li+p, Q = -4.02MeV.(4) 

For energies of projectiles ~10 MeV, the cross sections for these nonthermal processes 
are on the order of 100 mbn, and indeed 10 7 times larger than the SBBN cross section for 
producing 6 Li . This enhancement figure underlines the 6 Li sensitivity to non-thermal 
BBN, and makes it an important probe of energy injection mechanisms in the early 
Universe. 

Another important aspect of the nonthermal BBN is the possibility to alleviate the 
tension between the Spite plateau value and the predicted abundance of 7 Li, e.g. "solve 
the 7 Li problem". To achieve that the energy injection should occur in the temperature 
interval 60 keV > T > 30 keV, i.e. during or just after 7 Be synthesis . The essence of this 
mechanism consists in the injection of ^ 10~ 5 neutrons per baryon that will enhance 
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Figure 2. Upper bound on the annihilation cross section of particle dark 
matter of mass M x from BBN. Here the upper line assumes annihilation into 
only electromagnetically interacting particles, whereas the two lower lines assume 
annihilation into a light quark-anti-quark pair. Adopted limits on the light element 
abundances are as indicated in the figure. 



Be— > Li interconversion followed by the p-destruction of Li via the thermal reaction 
sequence [18] : 



Note that this is the same mechanism that depletes 7 Be in SBBN, but with elevated 
neutron concentration due to the hadronic energy injection. This mechanism of 
depleting 7 Be is tightly constrained by the deuterium abundance, as extra neutrons 
could easily overproduce D. 

Fig. [1] summarizes constraints on abundance vs lifetime of relic decaying particles. 
It is convenient to measure the abundance in terms of Vtxh 2 , the present day fraction of 
total energy density if these particles were to remain stable. This quantity relates to the 
in the literature frequently used ( = nxMx/s via ( = 3.6639 x 10~ 9 GeVQxh 2 where 
nx is particle number density, Mx is particle mass, and s is entropy. It is seen that 
constraints get increasingly more stringent when the lifetime Tx increases, implying also, 
that under generic circumstances, the production of dark matter X (with Vtxh 2 ~ 0.1) 
by the decay of a parent particle Y — > X + ... at r x 3> 10 3 sec is extremely problematic, 
if at all possible. 

5. Residual dark matter annihilation during BBN 

Many dark matter candidates X may be ultimately visible in our Galaxy due to the 
cosmic rays they inject induced by residual XX self-annihilations. This is, for example, 
the case for supersymmetric neutralinos, provided that their annihilation products 



n + 



Be -> p + 7 Li; p + 7 Li -> 4 He + 4 He . 



(5) 
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Figure 3. Dark matter annihilation rate versus dark matter mass. The blue band 
shows parameters where 6 Li due to residual dark matter annihilation may account 
for the 6 Li abundance as inferred in HD84937 ( 6 Li/ 7 Li« 0.014 - 0.09 at 2-a), 
whereas the orange-red-green-yellow region shows where 7 Li is efficiently destroyed 
i.e. 7 Li/H< 1.5,2,3, and 4 x 10~ 10 , respectively. Above the lower (upper) dashed 
line D/H exceeds 4 x 10~ 5 (5.3 x 10~ 5 ), such that parameter space above the upper 
dashed line is ruled out by D overproduction. Scenarios between this line and the 
upper edge of the blue band are problematic since severely overproducing 6 Li. Dark 
matter annihilation into light quarks has been assumed. 



can be distinguished from astrophysical backgrounds. Residual annihilation events 
in the early Universe may also be of importance as they may lead to cosmologically 
significant 6 Li abundances [52] induced by the non-thermal nuclear reaction discussed 
in Section HI Given an X annihilation rate (av) and X-density nx one may determine 
the approximate fraction fx of X particles which annihilate in the early Universe at 
temperature T 



where Atn ~ H^ 1 = (QOMpJ^gT^) 1 / 2 is the characteristic Hubble time at T, g is 
the appropriate particle statistical weight, s = 47r 2 /90gT 3 is radiation entropy, and 
(...) denotes a thermal average, which can be taken once the velocity dependence of 
av is specified. Many scenarios for the production of dark matter envision a stable 
self-annihilating particle, typically a weakly-interacting massive particle or WIMP, 
whose final asymptotic abundance is given by its annihilation rate. Straightforward 
considerations of thermal WIMP freeze-out require the annihilation rate at Tj h ~ 
0.05mx to be (avYJ 1 ~ lpbnxc = 3x 10~ 26 cm 3 s _1 if the X-particle is to be the dominant 
component of dark matter, VLxh 2 « 0.1. Less straightforward but still plausible scenarios 
that include the non-thermal production of dark matter, e.g. via evaporation of Q-balls 
and/or decay of relic particles Y — > X + ... with subsequent X self-annihilation, may 
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require (o~v)f well in excess of (avYJ 1 for Qxh 2 ~ 0.1. We concentrate on the WIMP 
example, and parametrize the velocity dependence of (av) by (av) = (av) S(v), with 
the chosen normalization (av)o = (o~v)f. We are interested in finding the fraction of 
annihilating WIMP particles at T < lOkeV, a temperature scale below which 6 Li is 
no longer susceptible to nuclear burning (destruction). Exploiting ([6]) at the freeze-out 
temperature Tf, where fx(Tf) ~ 1, as well as at an arbitrary other T <C Tf, we obtain 



for fx *C 1. Several generic options are possible for the temperature scaling of the 
the .S-ratio in (jTJ). If the s-wave annihilation is mediated by short-distance physics and 
occurs away from sharp narrow resonances, (S(v))t = (S{v))T f = 1, where the second 
equality is due to our chosen normalization. Using this conservative assumption and 
Eq. (ED, for a WIMP of mass M x = 100 GeV, so that g(T)/g(T f ) ~ 0.1, one finds 
that only a small fraction, fx ~ 6 x 10~ 7 , of X-particles has a chance to annihilate at 
T ~ lOkeV and below. Nevertheless, even this tiny fraction is still sufficient to produce 
a 6 Li abundance of 6 Li/H^ 1.6 x 10~ 12 . 

Existence of attractive Coulomb-like force of some strength a' in the WIMP sector 
may lead to a significant enhancement of annihilation at low temperatures/velocities [53J, 
possibly leading to a much higher yield of 6 Li. In this case the Sommerfeld-like scaling 
av ~ (ira'/v)[l — exp(— ira'/v)} -1 , enhances the annihilation at small v^na', i.e. 
(S(v))t — net' /v. This leads to a ~ T -1 / 2 scaling of (S(v))t in (Cj) when X-particles are 
still in thermal equilibrium with the plasma. After they have dropped out of thermal 
equilibrium, (S(v))t falls even more rapidly as ~ T" 1 , with the net effect that weak 
mass scale X-particles usually have much smaller velocities at the end of BBN than 
in the Milky Way. Similarly, the presence of narrow resonances just above the XX 
annihilation threshold may drastically boost the annihilation at low energies. Both 
mechanisms of enhancing the annihilation have been widely discussed (see e.g. |54|) 
in an attempt to link some cosmic-ray anomalies to dark matter annihilation, as for 
example an elevated positron fraction e + /(e~ + e + ) observed by PAMELA instrument 

Keeping the annihilation rate as a free parameter, Fig. [2] shows the upper limit on 
the effective annihilation cross section imposed by BBN. Here electomagnetically- (upper 
line) and hadronically- (lower lines) annihilating particle DM has been considered. The 
former is mostly constrained by overproduction of 3 He/D at T ~ 0.1 keV, while the 
latter by 6 Li/ 7 Li overproduction at T « 10 keV, such that the effective annihilation cross 
section refer to (erf) at those temperatures. Due to the possibility of 6 Li destruction 
a fairly conservative 6 Li/ 7 Li< 0.66 constraint has also been considered. It is seen 
that much 6 Li may be produced by hadronic annihilations. Fig. [3] shows dark matter 
parameters which lead to the production of a 6 Li/ 7 Li ratio as claimed to be observed 
in the star HD84937 at la and 2a (dark blue and light blue), respectively. Here a 
completely hadronic XX — > qq annihilation has been employed, an assumption which 
could be further confronted with the constraints on antiproton fluxes in our galaxy. 
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Electromagnetic annihilations, as often implied in recent dark-matter interpretations of 
PAMELA |55j. may also lead to 6 Li, but only annihilations below T<0.3keV may do 
so efficiently Note, that the figure already implicitly assumes a factor ~ 3 — 4 stellar 
destruction of 7 Li (and 6 Li) to solve the lithium problem. It is therefore found that 
weak-scale mass dark matter particles, if fairly light, and if annihilating into hadronically 
interacting particles, may account for all of the observed 6 Li in HD84937. The figure 
also shows by the orange-red-green-yellow areas the dark matter parameters which 
would lead to a significant 7 Li destruction due to residual dark matter annihilations, 
with the nuclear destruction mechanism described in the previous section. Since those 
regions are much above the 6 Li band, the possibility of a factor of 2 or more depletion 
in 7 Li is severely constrained by 6 Li overproduction. We note in passing that our 
constraints are significantly more conservative than those found in Ref. [56]. It is 
intriguing to realize that primordial 6 Li production by residual (hadronic) dark matter 
annihilations dominates standard BBN 6 Li production for annihilation rates as small as 
(av) ~ 10 _27 cm 3 s _1 , well below (avYJ 1 . It is thus possible that in the most primeval 
gas clouds and in the oldest stars the bulk of 6 Li is due to dark matter annihilations. 
Unfortunately, 6 Li abundances as low as 6 Li/H < 10~ 12 are difficult to observe. 

6. Catalyzed BBN (CBBN) 

The idea of particle physics catalysis of nuclear reactions goes back to the 1950s, and 
muon- catalyzed fusion has been a subject of active theoretical and experimental research 
in nuclear physics. In recent years there has been a significant interest towards a 
possibility of nuclear catalysis by hypothetical negatively charged particles that live 
long enough to participate in nuclear reactions at the BBN time [53, EH EH1 EH EH EH 
E21 E31 Ei ESI ESI EHl E7j (see also Ref. [6HI E3 [70] for earlier work on the subject). An 
essence of the idea is very simple: a negatively charged massive particle that we call 
X~ gets into a bound state with the nucleus of mass and charge Z, forming a large 
compound nucleus with the charge Z — 1, mass Mx + m^, and the binding energy in 
the 0(0.1 — 1) MeV range. Once the bound state is formed, the Coulomb barrier is 
reduced signalling a higher "reactivity" of the compound nucleus with other nuclei. But 
what proves to be the most important effects of catalysis, are new reaction channels 
which may open up and avoid SBBN-suppressed production mechanisms [57], e.g. Eq. 
(J21), thus clearing path to synthesis of elements such as 6 Li and 9 Be. Although in this 
chapter we discuss the catalysis by negatively charged heavy relics, this is not the only 
option for CBBN, as for example, strongly interacting relics may also participate and 
catalyze certain nuclear reactions. 

Although the connection between dark matter and CBBN is not immediate - after 
all the dark matter may not be charged - it is possible that dark matter particles 
do have a relatively long-lived charged counterpart. One example of this kind is 
supersymmetry with the lightest supersymmetric particle (LSP) the gravitino and the 
next-to-LSP (NLSP) a charged slepton to be examined in the next section. In that case 
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bound state 


a [tm) 


\E b \[keV] 


T [keV] 


P x- 


29 


25 


0.6 


4 HeX" 


3.63 


346 


8.2 


7 BeX" 


1.03 


1350 


32 


8 BeX- 


0.91 


1430 


34 



Table 1. Properties of the bound states: Bohr radius ao ~ l/(ZarriN), binding 
energies Et, calculated for realistic charge radii, and "photo-dissociation decoupling" 
temperatures To. 

the decay of the NLSP is tremendously delayed by the smallness of the gravitino-lepton- 

i li 

slepton coupling ~ . Another example in the same vein is the nearly degenerate 
stau-neutralino system, in which case the longevity of the charged stau against the 
decay to the dark matter neutralino is ensured as long as the mass splitting of the 
stau-neutralino system is below 100 MeV. Both, the gravitino and neutralino in these 
two examples represent viable dark matter candidates. A very important aspect of 
CBBN is that the abundance of charged particles before they start decaying is given by 
their annihilation rate at freeze-out. In most of the models their abundance is easily 
calculable, and if no special mechanisms are introduced to boost the annihilation rate, 
the abundance of charged particles per nucleon is not small, and in the typical ballpark 
of Y x ~ (0.001 - 0.1) x m x /TeV. 
Properties of the bound states. 

For light nuclei participating in BBN, we can assume that the reduced mass of the 
nucleus- X~ system is well approximated by the nuclear mass, so that the binding energy 
is given by Z 2 a 2 m J 4/2 when the Bohr orbit is larger than the nuclear radius. It turns 
out that this is a poor approximation for all nuclei heavier than A = 4, and the effect 
of the finite nuclear charge radius has to be taken into account. In Tabled] we give the 
binding energies, as well as the recombination temperature, defined as the temperature 
at which the photodissociation rate of bound states becomes smaller than the Hubble 
expansion rate. Below these temperatures bound states are practically stable, and the 
most important benchmark temperatures for the CBBN are the T ~ 30, 8, 0.5 keV, 
when ( 7 BeX~), ( 4 HeX~), and (joX~) can be formed without efficient suppression by the 
photodissociation processes. It is important to emphasize that these properties of the 
bound states are generic to any CBBN realization: i. e. they are completely determined 
by the charge of X~ and electromagnetic properties of nuclei, and thus are applicable 
to SUSY or non-SUSY models alike. It is also important to note that the ( 8 BeX~) 
compound nucleus is stable, which may open the path to synthesis of A > 8 elements 
in CBBN. 

Catalysis at 30 keV: suppression of 7 Be. 

When the Universe cools to temperatures of 30 keV, the abundances of deuterium, 
3 He, 4 He, 7 Be and 7 Li are already close to their freeze-out values, although several 
nuclear processes remain faster than the Hubble rate. At these temperatures, a 
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negatively charged relic can get into bound states with 7 Be and form a ( 7 BeX~) 
composite object. Once this object is formed, some new destruction mechanisms for 
7 Be appear. For models with weak currents connecting nearly mass-degenerate X~-X° 
states, a very fast internal conversion followed by the p-destruction of 7 Li: 

( 7 BeX~) -> 7 Li + X°; 7 Li + p -> 2a. (8) 

When X~ —>X is energetically disallowed the destruction of 7 Be can be achieved via 
the following chain: 

( 7 BeX") + p -> ( 8 BX~) + 7 : ( 8 BX~) -> ( 8 BeX~) + e + z/, (9) 

which is much enhanced by the atomic resonances in the ( 7 BeX~) system [SI] . 

The rates for both mechanisms may be faster than the Hubble rate, possibly 
leading to a sizable suppression of 7 Be abundance z/( 7 BeX _ ) bound states are efficiently 
forming. In other words, ( 7 BeX~) serves as a bottleneck for the CBBN depletion of 7 Be. 
The recombination rate per 7 Be nucleus leading to ( 7 BeX~) is given by the product of 
recombination cross section and the concentration of X~particles. It can be easily shown 
that for Y x < 0.01 the recombination rate is too slow to lead to a significant depletion of 
7 Be. Detailed calculation of recombination rate and numerical analyses of the CBBN at 
30 keV [6U [66] find that the suppression of 7 Be by a factor of 2 is possible for Yx > 0.1 
if only mechanism ([9]) is operative, and for Yx > 0.02 if the internal conversion (jSJ) is 
allowed. 

Catalysis at 8 keV: enhancement o/ 6 Li and 9 Be. 

As the Universe continues to cool below 10 keV, an efficient formation of ( 4 HeX~) 
bound states becomes possible. With the reasonable assumption of Y x < F H e the rate 
of formation of bound states per X~ particle is given by the recombination cross section 
and the concentration of the helium nuclei. Numerical analysis of recombination reveals 
that at T ~ 5 keV about 50% of available X~ particles will be in bound states with 
4 He[EZ|. 

As soon as ( 4 HeX~) is formed, new reaction channels open up. In particular, a 
photonless thermal production of 6 Li becomes possible 

( 4 HeX~) +D -f 6 Li + X~; Q ~ 1.13MeV, (10) 

which exceeds the SBBN production rate by ~six orders of magnitude. The production 
of 9 Be may also be catalyzed, possibly by many orders of magnitude relative to the 
SBBN values, with the following thermal nuclear chain [65] : 

( 4 HeX") + 4 He — > ( 8 BeX") + 7; ( 8 BeX~) + n -> 9 Be + X~ . (11) 

Both reactions at these energies are dominated by the resonant contributions, although 
the efficiency of the second process in (ITT]) is not fully understood. 

Current estimates/calculations of the CBBN rates are used to determine the generic 
constraints on lifetimes/abundances of charged particles. The essence of these limits is 
displayed in Figure HJ which shows that for typical X~ abundances the lifetime of the 
charged particles would have to be limited by a few thousand seconds! This is the main 
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conclusion to be learned from CBBN. Note that while, to lowest order, non-thermal 
BBN is sensitive to the energy density of decaying particles, the CBBN processes are 
controlled by the number density of X~, which underlines the complimentary character 
of these constraints. In some models, where both catalysis and cascade nucleosynthesis 
occur, catalysis dominates cascade production of 6 Li for all particles with hadronic 
branching ratio B h < 1CT 2 [51] . whereas 7 Li destruction is usually dominated by cascade 
effects unless Bh ^ 10~ 4 . 9 Be production, on the other hand, is conceivable only through 
catalysis. 

Catalysis below 1 keV and nuclear uncertainties 

Finally we comment on the possibility of {pX~) catalysis of nuclear reactions, 
discussed in Refs. [691 E3] • Although it is conceivable that the absence of the Coulomb 
barrier for this compound nucleus may lead to significant changes of SBBN/CBBN 
predictions, in practice it turns out that in most cases (pX~)-related mechanisms are 
of secondary importance. The large radius and shallow binding of this system leads to 
a fast charge-exchange reaction on helium, {pX~) + 4 He— > ( 4 HeX~) + p, that reduces 
the abundance of (pX~) below 10~ 6 relative to hydrogen, as long as Y x - "Z Y 4 He, making 
further reactions inconsequential for any observable element [38J. In the less likely case, 
Y x - ^^4 He , significant late-time processing due to (pX~) bound states may still occur. 
Such late time BBN, nevertheless, typically leads to observationally unacceptable final 
BBN yields. 

Unlike in the SBBN case and even in cascade nucleosynthesis that utilizes mostly 
measured nuclear reaction rates, CBBN rates cannot be measured in the laboratory, 
and significant nuclear theory input for the calculation of the reaction rates is 
required. However, since the X~ participates only in electromagnetic interactions, such 
calculations are feasible, and dedicated nuclear theory studies [60] in this direction has 
already commenced. The reaction rates for some CBBN processes, such as and (|T0l) 
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Figure 5. Parameter space in the relic decaying particle abundance times hadronic 
branching ratio Bh, i.e. Qx^Bh, and life time tx plane, where 7 Li is significantly 
reduced (red and blue) and 6 Li is efficiently produced (green and blue). See text for 
further details. From Ref. [7Tj . 



are already known within a factor of 2 accuracy, and the detailed calculations for the 
9 Be synthesis are underway 



7. Dark Matter Production during BBN: NLSP^LSP example 

Dark matter particles may be produced by the decay of relic parent particles X during 
BBN. Examples, well-studied by different groups, include the production of gravitino- 
LSP dark matter by NLSP decays (often charged sleptons or neutralinos) or production 
of neutralino dark matter by heavier gravitinos. Other conceivable possibilities include 
the production of superweakly interacting Kaluza-Klein dark matter, and more generally 
the cascade decays to any superweakly interacting dark matter candidates. In case of 
charged NLSP decays, both nonthermal and CBBN processes must be accounted for. 
In the framework of gravitino-LSP/stau-NLSP the lifetime of the charged slepton in the 
limit of rriQ «C M NL sp is given by 



TNLSP 



24xl0 4 sec x ( MnLSP V V ^ V 
2.4 xlU sec x ^ 300Gey j ^ 10GeV J 



(12) 



where M NLSP and ttIq denote NLSP and gravitino mass, respectively. 

It becomes exceedingly more difficult with increasing tx to obtain observational 
consistency with inferred primordial abundances (cf. Fig. [T|). BBN therefore plays an 
important role in constraining such scenarios (cf. Sections H] and E]) • However, BBN may 
not only constrain, but also favor particular scenarios, if current discrepancies with 6 Li 
and 7 Li abundances are to be taken seriously. Both trends, the reduction of 7 Li and the 
production of 6 Li via mechanisms described in Sections H] and El are seen in Fig. [5j There 
the red area shows decaying particle parameter space resulting in more than a factor of 
2 suppressed 7 Li abundance relative to the SBBN prediction, and the green area shows 
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Figure 6. Parameter space in the CMSSM which may impact the primordial 7 Li 
and/or 6 Li abundances with the color coding essentially as in Fig. [5] Left: unifying 
scalar soft mass Too versus gaugino soft mass TO1/2; Right: gravitino abundance fl^h 2 
versus NLSP life time r. ^From Ref. [7T] . 

regions where significant 6 Li production (0.015 < 6 Li/ 7 Li ^0.3) occurs. In the overlap 
of these areas, the blue area, both effects may be achieved simultaneously [1H1 EU 172] . 
Fig. [5] also shows the prediction of supersymmetric scenarios with the gravitino- LSP, for 
some representative values of other supersymmetric mass parameters. In particular, the 
grey dots show predictions of stau NLSPs with gravitino LSPs of mass tUq = 50 GeV 
within the so-called constrained minimal supersymmetric SM (CMSSM), whereas the 
blue dots show the case of neutralino NLSPs decaying into tjiq = 100 MeV gravitino 
LSPs within the gauge-mediated supersymmetry breaking scenario. It is seen that both 
scenarios naturally cross the region of " 7 Li destruction". The assumption underlying 
these models is a thermal freeze-out abundance of the NLSP. Since this typically leads 
to NLSP abundances, 10 -3 < I^nlsp ~ 10 3 , and taking into account that gravitino energy 
density due to NLSP decays is Qq = f2 NL sp (^g/^nlsp); the resulting Qq produced in 
such scenarios may come close to the observed dark matter density. This is particularly 
the case for heavy gravitinos itlq ~ 100 GeV in the CMSSM, for which a more detailed 
results are shown in Fig. [6j It is intriguing, and perhaps purely coincidental, that when 
resolving the tension between observed and predicted 7 Li abundance by staus decaying 
into gravitinos, the resulting gravitino abundance may account for all the dark matter. 
For stau decay times r pa 10 3 sec it is furthermore possible to synthesize a primordial 
6 Li abundance as claimed to be observed in low-metallicity stars. Moreover, though less 
certain, the same parameter space could also lead to an important 9 Be abundance due 
to catalytic effects (cf. Section [6]), as indicated by the cross-hatched region in Fig. El 
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Finally, since produced by decays, gravitino dark matter in such scenarios is significantly 
warm, with free-streaming velocities of the same order as those of a m ~ 3keV early 
freezing-out relic particle, which has important implications for the small scale structures 
in the present day Universe. It is therefore not impossible that some time in the future, 
anomalies in the primordial light elements may have been understood as signs of the 
dark matter. Nevertheless, independent verification by particle accelerators, such as the 
LHC is required. Unfortunately, scenarios as presented in Fig. [6] require staus of mass 
mf > 1 TeV too heavy to be produced at the LHC. 

8. Conclusions 

Even though the concept of Big Bang Nucleosynthesis is more than 50 yr old, it continues 
being relevant due to the constant progress in nuclear physics, astrophysics, and the 
refined quality of calculations. In this paper, we have reviewed the status of BBN, and 
have shown how New Physics can modify the synthesis of light element abundances. All 
three generic ways, extra degrees of freedom modifying the Hubble expansion during 
BBN, energy injection due to annihilation or decay of heavy particles, and particle 
catalysis of BBN reactions, are highly relevant to the physics associated with particle 
dark matter, or with particles intimately tied to dark matter. We have illustrated 
how the existing overall concordance between the predicted elemental abundances and 
observations lead to some very non-trivial constraints on the properties of DM particles 
and their companions {e.g. stau-gravitino system). 

Perhaps even more intriguing is the current discrepancy between the observed and 
standard BBN predicted abundance of 7 Li at the level of 2-3. This discrepancy has 
firmed up since the last unknown SBBN parameter, the baryon-to-photon ratio, has been 
determined with better than 5% acuracy by recent high-precision CMB experiments. At 
this moment it is premature to tell how the 7 Li problem is resolved, but it is nonetheless 
intriguing that certain models with unstable particles are capable of alleviating this 
discrepancy. Hopefully, the continuing improvement of observational determination of 
primordial light element abundances, as well as the future breakthroughs in electroweak 
scale particle physics would help to solve this important problem. 
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